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1. Introduction 

The QCD interaction is flavour-blind. Neglecting electromagnetic and weak interactions, the 
only difference between quark flavours comes from the quark mass matrix. We investigate here how 
flavour-blindness constrains hadron masses after flavour SU (3) is broken by the mass difference 
between the strange and light quarks, to help us extrapolate 2 + 1 flavour lattice data to the physical 
point. 

We have our best theoretical understanding when all 3 quark flavours have the same masses 
(because we can use the full power of flavour SU (3)); nature presents us with just one instance of 
the theory, with ml / ml k, 25. We are interested in interpolating between these two cases. 

We consider possible behaviours near the symmetric point, and find that flavour blindness is 
particularly helpful if we approach the physical point, (m^ * , * ) , along a path in the {m^ , mf ) plane 
starting at a point on the SU (3) flavour symmetric line (mf = mf) and holding the sum of the quark 
masses mf, + m^, + = 2m^ + constant, as sketched in Fig. 




Figure 1: Sketch of the path (red, solid line) in the {mf,m'^) plane to the physical point {mf*,m^*). 



2. Theory 

Our strategy is to start from a point with all 3 sea quark masses equal, 

m'^ = m''^=m^=m'o, (2.1) 

and extrapolate towards the physical point, keeping the average sea quark mass 

m' = ^,{m: + m''^ + m'^) (2.2) 

constant. For this trajectory to reach the physical point we have to start at a point where « I'wf *. 
As we approach the physical point, the u and d quarks become lighter, but the s becomes heavier. 
Pions are decreasing in mass, but K and rj increase in mass as we approach the physical point. 
We introduce the notation 

5m'^ = m^^—Wf , q = u,d,s, (2.3) 
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and later use a similar notation for bare quark masses. (We will be mainly interested in the 2 + 1 
flavour case, with m"=m'^ = ml.) With this notation, the quark mass matrix is 



m 



m 






(l 


o\ 
















-1 







-1 








J 












(2.4) 



The mass matrix ^ has a singlet part (proportional to /) and an octet part, proportional to A3, 
Ag. In the 2+1 case = Sni^ and isospin is a good symmetry. We argue that the theoretically 
cleanest way to approach the physical point is to keep the singlet part of constant, and vary 
only the non-singlet parts. One technical advantage of this strategy is that it simplifies the quark 
mass renormalisation. In the case of cloverAVilson fermions, the singlet and non-singlet parts of 
the mass matrix will renormahse with different renormalisation constants [Q] 



Z^^img + azm) 



(2.5) 



where az represents the fractional difference between the renormalisation constants. (Numerically 
this factor is 0(1), and is thus non-negligible. Of course, for chiral fermions az = 0.) This gives 



m 



:Z'Jil + az)m, 



(2.6) 



and so by keeping the singlet mass constant we avoid the need to use two different Zs. This means 
that even for clover actions it does not matter whether we keep the bare or renormalised average 
sea quark mass constant (so we shall drop the " superscript in the following considerations). 

An important advantage of our strategy is that it strongly constrains the possible mass de- 
pendence of physical quantities, and so simplifies the extrapolation towards the physical point. 
Consider a flavour singlet quantity (for example the scale ro, or the plaquette P) at the symmetric 
point {mo, mo, mo). If we make small changes in the quark masses, symmetry requires 



dro dro dvQ 



1m„ 



dmd 



im. 



If we keep m„ + + constant, dm^ = —dmu — dm^j = —2dmi so 



dro 



dro dro , , ^''o 
: dm,, h dmij h dm,,. - — 



m,, 



dm. 



im.. 



0. 



(2.7) 



(2.8) 



The effect of making the strange quark heavier exactly cancels the effect of making the light quarks 
lighter, so we know that ro must have be stationary at the symmetrical point. This makes extrapo- 
lations towards the physical point much easier, especially since we find that in practice quadratic 
terms in the quark mass expansion are very small. Any permutation of the quarks, such as an inter- 
change M o or a cyclic permutation m — )• J — )• s — )• m doesn't change the physics, it just renames 
the quarks. Any quantity unchanged by all permutations will be flat at the symmetric point, like 
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Figure 2: The behaviour of the octet and decuplet under the permutation group Sj,. The colours denote sets 
of particles which are invariant under permutations of the quark flavours. 



ro- We can also construct permutation-symmetric combinations of hadrons. For example, for the 



decuplet, any permutation of the quark labels will leave the £ , (uds) unchanged, so the £ is 
shown by a single black point in Fig. ^ On the other hand, a permutation (such as u ^ d ^ s) 
can change a A++(mmm) into a A^{ddd) or (if repeated) into an D.^{sss), so these three particles 
form a set of baryons which is closed under quark permutations, and are all given the same colour 
(red) in Fig. ^. Finally the 6 baryons consisting of two quarks of one flavour, and one quark of a 
different flavour, form an invariant set, shown in blue in Fig. ^. If we sum the masses in any of 
these sets, we get a flavour-symmetric quantity, which will obey the same argument we gave in 



eq. ( [2.8| ) for the quark mass (in)dependence of the scale tq. We therefore expect that the £ mass 



must be flat at the symmetric point, and furthermore that the combinations (Ma++ +Ma- +Ma) 
and (A/a+ + M^o + Mj;*+ + M^*- + M;^*o + M^.- ) will also be flat. TechnicaUy these symmetrical 
combinations are in the A[ singlet representation of the permutation group ^ ^3. We list some of 
these invariant mass combinations in Table |l|. We can use the singlet combinations from this table 



Decuplet 


2MA + Ma 


red 


baryons 


2(Ma + Mz*+Ms.) 


blue 






black 


Octet 


2{Mn + M^+Mz) 


blue 


baryons 


Mi;+Ma 


black 


Pseudoscalar 


4MI + 2MI 


blue 


mesons 




black 


Vector 


AMK'+2Mp 


blue 


mesons 


2Mp+M^^ 


black 



Table 1: Permutation invariant mass combinations, see Fig. g. (j)s is a fictitious ss particle; Tjg a pure octet 
meson. The colours in the third column correspond to Fig. |[ 



to locate the starting point of our path to physics by fixing a dimensionless ratio such as 



4- = physical value , 



(2.9) 



N 



'53 has the same symmetry group as that of an equilateral triangle, C3,.. This group has 3 irreducible representations, 
pp, two different singlets, A; and A2 and a doublet E, with elements £+, . 
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where = ^{IM^^+Ml) and = ^{Mn + Mz + Mz). The permutation group ^3 yields a lot of 
useful relationships, but cannot capture the entire structure. For example, there is no way to make 
a connection between the A++(mmm) and the A^{uud) by permuting quarks. To go further, we need 
to classify physical quantities by SU (3) and the permutation group ^3 (which is a subgroup of 
SU{3)). 

Let us first consider hnear terms in Snig. These are given in Table |2[ Since we are keeping 



Polynomial 




S3 


SU{3) 
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Ai 
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(m — mo) 




Ai 
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5ms 
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E+ 
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{dniu - dnid) 


/ 


E 




8 


(rn — mo)^ 




Ai 


1 




{m — mo)5ms 




E+ 




8 


(m-mo)(5m„-5mrf) 




E 




8 


5ml + 5m^ + 5m^ 


/ 


Ai 


1 


27 


35m^ — {8mu — 8miiY 


/ 


E+ 




8 27 


5ms{dmd - dmu) 


/ 


E 




8 27 



Table 2: All the quark-mass polynomials up to (9(m^), classified by symmetry properties. A tick (/) marks 
the polynomials relevant on a constant m surface. These polynomials are plotted in Fig. ^. 



m constant, we are only changing the octet part of the mass matrix in eq. (2.4). Therefore, to 



first order in the mass change, only octet quantities can be effected. SU (3) singlets have no linear 



dependence on the quark mass, as we have already seen by the symmetry argument eq. ( |2.8| ), but 
we now see that all quantities in SU (3) multiplets higher than the octet cannot have linear terms. 
This provides a constraint on the hadron masses within a multiplet and leads to the Gell-Mann 
Okubo mass relations [Q]. 

In the 2 + 1 limit the decuplet baryons have 4 different masses (for the A, £*, E*, and H), but 
there is only one slope parameter in the Unear mass formula. Similarly, for the octet baryons there 
are 4 distinct masses, (A'^,A,r,E), but only 2 slopes; and for octet mesons, one slope parameter for 
three mesons, {n,K,ri). Mesons have fewer slope parameters than baryons because of constraints 
due to charge conjugation. In the meson octet the K and K must have the same mass, but there is 
no reason why the N and S, (which occupy the corresponding places in the baryon octet), should 
have equal masses once flavour SU (3) is broken. 

When we proceed to quadratic polynomials we can construct polynomials which transform 
like mixtures of the 1, 8 and 27 multiplets of SU (3), Table ^ This covers all the structures that can 
arise in the octet mass matrix, but the decuplet mass matrix can include terms with the symmetries 
10, 10, and 64, which first occur when we look at cubic polynomials in the quark masses. Table ^. 
The allowed quark mass region on the m = const, surface is an equilateral triangle, as shown in 
Fig. ^. Plotting the polynomials of Tables ^, ^ then gives the figures in Fig. ^ where the colour 
coding indicates whether the polynomial is positive (red) or negative (blue). 
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Polynomial 
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E+ 




i 27 




[rn — mo)5ms{dmij — dm^) 




E- 
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i 27 




5mu5mii5ms 
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Ai 
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27 


64 


dms{dml + dm^ + Sm^) 
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E+ 




i 27 


64 


(5m„ - dmd) {5ml + Sm^ + dm^) 


/ 


E- 




i 27 


64 


{5ms - 5mu){5ms - 5md){5mu - dm^) 


/ 


Ai 




10 To 


64 



Table 3: The cubic quark-mass polynomials, classified by symmetry properties. A tick (/) marks the 
polynomials relevant on a constant Tn surface. These polynomials are plotted in Fig. Q 



(5 m 



(5 m 




(5m„ 



(5 m, 



Figure 3: The allowed quark mass region on the m = const, surface is an equilateral triangle. The black 
point at the center is the symmetric point, the red star is the physical point. 2+1 simulations lie on the 
vertical symmetry axis. The physical point is slightly off the 2+1 axis because m^/ > ot„. 



We can see how well this works in practice by looking at the physical masses of the decuplet 
baryons, 

4Ma + SMz* + 2Mz* +Ma= 13.821 GeV singlet oc (Sm/)" 

-2Ma +Ms*+M£i = 0.742 GeV octet oc (5m;)^ 

4Ma - SMz" - IMz" + 3Ma = -0.044 GeV 27plet oc {Sm,)^ 
-Ma + 3Mi:* - 3Me* +Ma = -0.006 GeV 64plet oc (Sm,)^ (2. 10) 

When we form combinations with particular SU (3) symmetries we see a strong hierarchy, which 
suggests a short Taylor series may work well all the way from symmetry point (mo, mo, mo) to the 
physical point. This gives the constrained fit formulae 

Ml = M^ + 2a5m, + (j8o + 2j8i)5m? 
mI = mI- aSmi + (jSo + 5jSi + 9l52)5mf 

M^-^ = M^-4a5m, + {po + Spi)dmj , (2.11) 
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Ai E+ E- Ai 




Figure 4: Contour plots of the polynomials relevant for the constant m Taylor expansion, see Tables g, ^. A 
red(dish) colour denotes a positive number while a blue(ish) colour indicates a negative number If Wj, = m^ 
(the 2+1 case), only the polynomials in the Ai and columns contribute. 

Mp = Mo + 2a5m, + (jSo + 2Pi) 5mj 
Mk' = Mo - admi + (j8o + SjSi + 9^2)8m^ 

M^, = Mo-4a5mi + {po + ^Pi)Smj, (2.12) 



Ma^ = Mo + 3Aidmi + (Bo + 3Bi)dmj 

Ma = Mo + dAjSmi + (Bo + 6Bi - 3^2 + 9B^)5m] 

Mz = Mo- 3A25mi + (Bo + 6Bi + 3B2 + 9B3)5m} 

Ms = Mo-3(Ai -A2)5m/ + (Bo + 9Bi -3B2 + 9B3)5m^, (2.13) 



Ma = Mo + 3A5mi + (Bo + 3Bi)dmj 
Mz' = Mo + + (Bo + 6Bi + 9B2)dmj 
Me' = Mo- 3A5mi + (Bo + 9Bi + 9B2)dmj 

Ma= Mo- 6A5m, + (Bo + \2Bi)5mj . (2.14) 



While the Unear terms are highly constrained, the quadratic terms much less so; indeed only for 
the baryon decuplet is there any constraint^ . Note also that for the pseudoscalar octet, M^,^ is a 
fictitious 55 particle (due to non-perfect tj-tj' mixing), while for the vector octet due to near perfect 



^The coefficients of the dmj term appear complicated; indeed there seem to be too many for the nucleon octet, 



eq ( 2.13 ). However although not discussed here the SU{3) flavour symmetry breaking expansion can be extended to 
different valence quarks than sea quarks or 'partially quenching'. In this case the m/, sea quarks remain constrained 
by m = const., but the valence quarks denoted by /i/, jls are unconstrained. With = jlq — m, the nucleon octet, 
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mixing between the ^ and ft) the is (almost) a perfect ss state, so that M^^ f« . However even 
these fictitious particles still contain useful information, due to the constrained fit. 



As eqs. ( [2.1 1[ )-( [2.14[ ) have been derived using only group theoretic arguments, then they will 
be valid for results derived on any lattice volume (though the coefficients are still functions of the 
volume). 

What is the effect of higher order terms on our path in the {ms,mi) plane? While at lowest 
order it did not matter whether we kept the quark mass constant or a particle mass now it does. 
Practically it is easiest to keep the (bare) quark mass fixed. Thus lines of constant Xs are curved 
though they do still have to have the slope of —2 at the point where they cross the SU (3) flavour 
symmetric line. Now we have to specify more closely what we mean when we keep 2m^ + 
constant, as different scale choices give slightly different paths as sketched in Fig. H. Note that the 



(m,,mj 



m =consJt-' 




Figure 5: Sketch of some possible paths (red lines) in the (ot/,OTj) plane to the physical point {mf*,m^*). 



physical domain m^,mj > translates to 



mi > 



■OCz 



(i + iotz 



-m. 



(l + Iaz) 



mi, 



(2.15) 



leading to non-orthogonal axes and possibly negative bare quark mass. (These features disappear 
for chiral fermions when az = 0.) If we make different choices of the quantity we keep constant at 
the experimentally measured physical value, for example 



y2 y2 



Y2 

xr 



(2.16) 



eq. (2.13) becomes 



Mn = M^) + 3AiS^li+BoSmj + 3BiS^lf 
Ma = Mi)+Ai{2S^ii + 8^,)-A2{5^,-8^,) 
Mz = Mf,+A^^{2S^l + 8^l,)+A2{8il,-8il,) 
Ms = M()+Ai{28^li + 8p.s)-A2{8ils-8ili) 



-BQ8mf+Bi{2S^f 
-Bo8mj+Bi{2Silf 
-Bo8mj+Bi{8^f- 



-8p.^)-B2{8ii^-8iif) 



-28tij) 



-B2{8tij-S^f) 



-B4{8^i,-8iiif 
-Bi{8fi,-8fiif 



-B2i8fij 



8ii}) 



-B3{8fi, 



-8ii,f. 



When fl ^ m (i.e. r eturn t o the 'unita ry lin e') then these results collapse to the previous results of eq. ( 2.13 ). Similar 
results hold for eqs. ( 2.11 ), ( [2.12| ) and ( 2.14 ). We can also show that on this trajectory the errors of the partially quenched 
approximation are much smaller than on other trajectories. 
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where in addition to the previously defined singlet quantities we also now have Xa = ^ (IM/i^+Mci), 
Xp = ^{2Mk' +Mp), = ^ we get slightly different trajectories. The different trajectories begin 
at slightly different points along the flavour SU (3) symmetric line. Initially they are all parallel 
with slope —2, but away from the symmetry line they can curve, and will all meet at the physical 
point. (Numerically we shall later see that this seems to be a small effect.) 

Finally we mention the connection of this method to ;tPT. For example for the pseudoscalar 
octet, using the results in ^ and assuming their validity up to the point on the flavour symmetric 
line, we find 



M, 



1 



^(3L4 + L5- 
/o 



-6L6-2L8) + 



X 



a 
Pi 



Qo 



16Z 



(3L4 + 2L5-6L6"4L8) + 



247r2/2 
X 



1 ^ 
X 



Ql 



Ql 



32(L5-2L8) + 



/32 = f 
Jo 



16(L5-2L8 




(2.17) 



with Qo = BqZ^jj, X = 2Qo{\ + a^)m, where the Ls are appropriate Low Energy Constants or LECs. 
For clover fermions, as mentioned before we have to respect the fact that singlet and non-singlet 
quark masses renormalise differently which leads to a non-zero az- 

We first note that when expanding the ;t-PT about a point on the SU (3) flavour symmetry 
Une gives to leading order only one parameter - a. (This means, in particular, that flavour singlet 
combinations vanish to leading order, as discussed previously.) Secondly, while we fit to a and jSq, 
j8i and jSa, it will be difficult to determine the individual LECs. The best we can probably hope for 
are these combinations. 



3. The Lattice 

We use a clover action for 2 + 1 flavours with a single iterated mild stout smearing as described 
in Also given in this reference is a non-perturbative determination of the improvement coeffi- 
cient for the clover term, using the Schrodinger Functional method. All the results given here will 
be at j3 = 5.50, c^w = 2.65. HMC and RHMC were used for the 2, 1 fermion flavours respectively, 
[^, to generate the gauge configurations. The (bare) quark masses are defined as 

' (3.1) 



2, yKq Ksym,c , 

where vanishing of the quark mass along the SU (3) flavour symmetric line determines Ksym\c- We 
then keep m = const. = wiq which gives 

= _^^_2 ■ (3-2) 
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So once we decide on Ki this then determines K^. A series of runs along the SU (3) flavour line 
determines the initial point on this line: K^ym by looking when Xj^/Xg, S = N,A, p are equal to their 



physical values, see eqs. ( [2.9[ ), ( |2.1q ). This would also include r if we have previously determined 
the physical value of tq. 

In Fig. ^ we plot X^/X^ for various Xg (with S = N, A, p, r). Also shown are constant 



0.30 



1 — ^ -1 



0.25 



0.20 



0.15 



0.10 



0.05 



0.00 



0.25 



0.50 



0.75 



i ^- 



• S=N 

• S = A 

• S = p 

• S = r 



1.00 



1.25 



Figure 6: X^/X^ versus M^/X^ for S ~ N, A, p, r for fciy,], = 0.12090. The dashed vertical line represents 
the physical value, while the dotted line gives the SU{3) flavour symmetric point. Filled points are on 
32^ X 64 lattices while open points are on x 48 sized lattice. Dashed horizontal lines represent constant 
fits to either the 32^ x 64 or 24^ x 48 results. The physical values are denoted by stars. For ro the small 
star is the value obtained using ro — 0.5 fm while the large star uses ro — 0.467 fm (determined for rif — 2 
degenerate flavours [^). 

fits to the two volumes — 24^ x 48 and 32^ x 64. (The runs on 24^ x 48 lattices have 0(2000) 
trajectories, while the runs on 32^ x 64 lattices have 0(1500) trajectories.) As mentioned before 
we first simulate at various quark masses on the SU (3) flavour symmetric line (i.e. m^/X^ = 1 in 
Fig. |6| where K^ym = 0.12090) to determine Ksym. Note that simulations with a Tight' strange quark 
mass and heavy 'light' quark mass are possible - here the right most point. (In this inverted strange 
world we would expect p or A decays.) It is apparent that while some scales fluctuate others 
(in particular S = N and r) are very smooth. We take this as a sign that singlet quantities are very flat 
and the fluctuations are due to low statistics. Assuming that all ratios are constant, (i.e. higher order 



effects are small, see eq. ( [2.16| )) we would expect all the ratios to converge to their experimental 
values. This is complicated because there are small finite size effects present (this can again best be 
seen in the S = N and r data). To investigate possible effects we are generating results on a variety 
of lattices, but for the present for the ratios we simply take the largest volume available. In Fig. ^ 
we plot aXs/Xs, S = N, A, p, r, n for the largest volume fitted results from Fig. ^ (together with 
smaller data sets for K^ym = 0.12080, K^ym = 0.12092, the latter data set is presently only partially 
generated). This ratio gives estimates for the lattice spacing a for the various scales. There appears 
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0.110 
0.105 
0.100 h 



— 0.095 
E 



5< 0.090 - 

X 

CD 

0.085 
0.080 h 
0.075 
0.070 



=0.12095 
lS = N 
>S = A 
I S = p 

I S = r [r„=0.467fm] 

lS = 7t 



8.268 



1 



8.270 



8.272 



8.274 



8.276 



8.278 



Figure 7: aXs/Xs against l/jciyi, for S = N, A, p, r, n and Kn™ 
dashed line is at K^yl = 0.12095. 



0.12080, 0.12090, 0.12092. The vertical 



to be convergence to a common scale of a ~ 0.078 fm. To investigate this point further we are 
performing additional runs at Ksym = 0.12092 and 0.12095. 



4. Spectrum results 

We now consider the mass spectrum. First we check whether there is a strong hierarchy due 
to the SU (3) flavour symmetry as found in eq. ( 2.10 ). In Fig. ^ we plot (4Ma + 3m^* + 2Mz* + 
Mn)/XA, {-2Ma + Ms* +M^)/Xa, (4Ma - SM^. - IMz- +Mci)/Xa and (-Ma + SM^. - SMs* + 
Mq)/Xa against M^/X|. Also shown are the experimental values. There is reasonable agreement 
with these numbers. (See [Q] for a similar investigation of octet baryons.) It is also seen that as 
expected while (— 2Ma +Ms* +Mq)/Xa has a linear gradient in the pion mass, in the other fits 
any gradient is negligible. To check for possible finite size effects we also plot a run at the same 
(k"/, Ks) but using 24^ x 48 lattice rather than 32^ x 64. There is little difference and so it appears 
that considering ratios of quantities within the same multiplet leads to (effective) cancellation of 
finite size effects. 

In Fig. ^ we plot the nucleon octet Mnq/Xn (for No = N, A, Z, S) against M'^/X^ for K^ym = 
0.12090. A typical 'fan' structure is seen with results radiating from the common point on the 
symmetric line. Again, finite volume effects tend to cancel in the ratio (normalising with the 
singlet quantity from the same octet) and so both volumes have been used in the fit. The combined 



fit uses eqs. (2.13), (2.11) with the bare quark mass being an 'internal' parameter. Note that one 
point has a light strange quark and a heavy 'light' quark. Similarly in Figs. |l^ and [Tl| we plot 
the corresponding nucleon decuplet Ma^/Xa and vector octet Mp^/Xp against M^/X^ respectively. 
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Figure 8: (4Ma + 3Mz* + 2Mz* +Ma)/XA, {-2Ma+Mz* +Ma)/X^, (4Ma - 5Mz* - 2Mz* +Mn)/XA 
and (-Ma + 3Mz* - 3Mh* +Mq) /X^ (filled circles) against mI/X^ together with a fit of constant, linear 
quadratic and cubic term in 5mi respectively. Extrapolated values are shown as opaque circles. Experimental 
values are denoted by stars. The opaque triangle is a run at the same (rc/, fCj), but on a 24^ x 48 lattice rather 
than a 32^ x 64 lattice. 
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Figure 9: Mn^/Xn (No = N, A, E, Z)) against mI/X^ together with the combined fit of eqs. ( |2.13| ), (gTTj) 
(the dashed lines). Experimental values are denoted by stars. The opaque triangle is a run at the same 
(k-/, K-.(), but on a 24^ x 48 lattice rather than a 32^ x 64 lattice. 
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Figure 10: Ma^/^a (Ad = A, E*, S*, £1) against M^/X^ together with the combined fit of eqs. ( [2.14| ), (|]T|) 
Same notation as Fig. 0. 
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Figure 11: Mp^/Xp (po — p, K*, 0s) against M^/X^ together with the combined fit of eqs. (2.12), (2.11) 
(the dashed lines). Same notation as Fig. 0. 
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Although we have included quadratic terms in the fit, there is really very little curvature in the 
results. We find good agreement with the experimental results. 

5. Conclusions 

We have outlined a programme to systematically approach the physical point starting from 
a point on the SU (3) flavour symmetric line. Exploratory results for the hadron mass spectrum 
show that constrained linear extrapolations give accurate results for the mass spectrum. We are 
also applying this method to the computation of matrix elements, some initial results are given in 



[10, 11] 



Acknowledgements 

The numerical calculations have been performed on the IBM BlueGeneL at EPCC (Edin- 
burgh, UK), the BlueGeneL and P at NIC (Juhch, Germany), the SGI ICE 8200 at HLRN (Berlin- 
Hannover, Germany) and the JSCC (Moscow, Russia). We thank all institutions. The BlueGene 



codes were optimised using Bagel, [|12|]. This work has been supported in part by the EU grants 
227431 (Hadron Physics2), 238353 (ITN STRONGnet) and by the DFG under contract SFB/TR 
55 (Hadron Physics from Lattice QCD). JZ is supported by STFC grant ST/F009658/1. 

References 

[1] W. Bietenholz et al, [QCDSF-UKQCD Collaboration], Phys. Lett. B690, 436 (2010) 

[arXiv: 1003.1114 [hep-lat] ]. 

[2] M. Gockeler et al, [QCDSF-UKQCD Collaboration], Phys. Lett. B639, 307 (2006) 

[arXiv:hep-ph/04 9312]. 

[3] P. W. Atkins, M. S. Child and C. S. G. Phillips, "Tables for Group Theory", (Oxford, 1970). 
[4] M. Gell-Mann, Phys. Rev. 125, 1067 (1962); S. Okubo, Prog. Theor Phys. 27, 949 (1962). 

[5] C. Allton et al. [RBC-UKQCD Collaboration], Phys. Rev D78, 1 14509 (2008) 

[arXiv: 0804 .04 73 [hep-lat] ]. 

[6] N. Cundy et al., [QCDSF-UKQCD Collaboration], Phys. Rev D79, 094507 (2009) 

[arXiv: 901 . 3302 [hep-lat] ]. 



[7] Y. Nakamura and H. Stuben, PoS (Lattice 2010)040 



[8] A. Ali Khan et al., [QCDSF Collaboration], Phys. Rev D74, 094508 (2006) 

[arXiv: hep-lat/ 60302 8]. 

[9] S. R. Beane et al., [NPLQCD Collaboration], Phys. Lett. B654, 20 (2007) 

[arXiv: hep-lat/ 604 013]. 



[10] J. Zanotti et al., [QCDSF-UKQCD Collaboration], [PoS (Lattice 2010) 165 
[11] F Winter etal, [QCDSF-UKQCD Collaboration] , |PoS (Lattice 2010) 16 3 
[12] P A. Boyle, Comp. Phys. Comm. 180, 2739 (2009). 



14 



